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Abstract: The Brodsky-Lepage-Mackenzie procedure is sequentially and unambiguously extended to 
any fixed order of perturbative QCD beyond the so called "large-/3o approximation" . As a result of this 
procedure, the obtained perturbation series looks like a continued-fraction representation. A subsequent 
generalization of this procedure is developed, in order to optimize the convergence of the final series, 
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Q_i' along the lines of the Fastest Convergence Prescription. This generalized BLM procedure is applied to 
O 1 the Adler D function and also to R e + e - in QCD at N 3 LO. A further extension of the sequential BLM is 
^ ■ presented which makes use of additional parameters to optimize the convergence of the power-series at 

r* ■ any fixed order of expansion. 
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1. Introduction 

The first goal of this article is to extend the well-known Brodsky, Lepage and Mackenzie (BLM) proce- 
dure of scale setting to any fixed order of perturbative QCD (pQCD). We will show that higher orders 
of pQCD in the MS-scheme unambiguously determine the new scales in the sequentially extended BLM 
prescription. More specifically, the effects of the coupling renormalization, encoded in the /^-function coef- 
ficients, are absorbed into a set of proper scales (if of the couplings = a s (/i?)/(47r) at any fixed order of 
pQCD. Thus, in contrast to the prevailing opinion, presented in the Particle Data Group booklet || (see 
rev. "Quantum Chromodynamics" therein), higher-order corrections can indeed "fix" the scales for the 
sequentially extended BLM (seBLM) procedure. At the same time, it may happen that this unambiguous 
seBLM choice of Oj, nf is not the "best" one from the point of view of the best series convergence, so that 
one should treat these two phenomena separately. For this reason, our first goal includes the analysis of 
the final seBLM result with an explicit demonstration of the above-mentioned problem. Our second goal 
following from this problem is to supply the seBLM procedure with a mechanism resembling the Fast 
Apparent Convergence (FAC) || (but not coinciding with it) in order to improve the convergence of the 
scries. 
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To simplify the analysis of the structure of the radiative corrections, renormalization- group invariant 
quantities, like the Adler function (D), are considered below. Quantities like the heavy quark potential Vq 
Hj], the Bjorken sum rules, the Gross-Llewellin Smith sum rule and so on, see, for example, Ref. |], || for 
a review, can be considered in the same manner. Here an empirical relation between the QCD /3-function 
coefficients bi, bi ~ b Q +1 has been used that will be explained later. The hierarchy of the contributions of 
the coupling renormalization to the perturbative coefficients at every order of a s = a s (fi 2 ) / '(47r) is based 
on this power law for b{. This detailed hierarchy requires a matrix representation for the perturbation- 
series expansion (PE) rather than the standard series one. The mentioned power relation works well 
at least up to the currently known coefficient 63 and Nf = 0-^5 active flavors in contrast to the 
usually discussed proposition of the so-called "large-60" limit that really implies — Nf ^> 1 ( therefore, 
bo ~ b\ in that approximation), see for a review Q. The seBLM procedure is based on this matrix 
representation; therefore, it can be formulated in terms of the dynamical characteristics, namely, the 
/3-function coefficients, rather than in terms of certain SU c (3)-Casimir operators that may appear at an 
intermediate stage. The use of the relation bi ~ 6q +1 allows one to go beyond this large-60 approximation 
in a natural way. Note that this power law should, of course, fail at some higher order of the PE when 
its expected factorial divergence starts to become operative. The proposed coupling renormalization is 
worked out (Sec. || and|6|) and includes as a particular case the "bubble approximation" elaborated on in 
|8|, |9|, The corresponding new perturbation series obtained with the seBLM procedure is also worked 
out. Then the /3-function expansion is performed explicitly for the 4-loop D-function and, subsequently, 
the seBLM procedure is applied to this D-function to highlight the advantages and disadvantages of the 
procedure in the case of this physical quantity. Although higher-order expansion coefficients have been 
calculated for other observables, mentioned above, the calculation of the beta-function expansion for 
these quantities constitutes a different task, and asks for a careful treatment. For this reason, we restrict 
ourselves here to the consideration of the D and R functions. 

Different kinds of extensions of the BLM approach have been discussed in a number of interesting 
articles [||, H [L(J 12, ||, 11, 13, which appeared in the last decade. The important issues about 



the scheme ambiguity of the BLM procedure [15], as well as the role of the anomalous dimensions (for 



the corresponding quantities) in the optimization procedures [16], will not be discussed here, with the 
corresponding calculations being performed in the MS-scheme for massless QCD. The present approach 
differs from those works mainly in three items: 

(i) it is based on a detailed matrix presentation for the PE; 

(ii) the proposed generalized scheme is formulated explicitly for any fixed order of the PE; 

(iii) all the sources of the coupling renormalization are taken into account and absorbed into the coupling 
scales. 

The seBLM procedure may not necessarily pertain to an improvement of the perturbation series. 
A relevant example was mentioned even in the pioneering work of Ref. An improved machinery, 
which uses the proper scales of the seBLM, applying a FAC-inspired optimization, is considered in Sec. ||. 
This generalization of the seBLM optimizes the perturbation series and formally completes the developed 
method in the sense of convergence. Because this method employs the apportions of the proper scales of 
the seBLM procedure in order to achieve the best convergence, it is named here xBLM procedure. 

This work is outlined as follows. In the next section, a convolution representation is proposed, which 
can be useful for the interpretation of the BLM procedure. This interpretation naturally leads to a 
certain kind of the "skeleton expansion" and is presented in a complete form in Appendix [C]. Appendix 
|A| , referred to in the same section, contains the proof of the convolution representation and the necessary 
formulas for the /3-function coefficients. Section ||] describes in detail the standard BLM procedure and 
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defines the task, while Sec. || presents the matrix rearrangement of the perturbation-series expansion. The 
first stage of the sequential generalization to a fixed iV-order case is considered in Sec. [|. The remaining 
steps and the completion of the task are given in Sec. ^, with some technical details presented for the 
convenience of the reader in App. [B]. An application of the developed method to the Adler function is 
considered in Sec. |7], with relevant technical details relegated to App. |D[ Finally, Sec. |8| presents the 
optimized version of our approach that yields a better convergence of the perturbation-series expansion. 



2. Convolution representation for amplitudes 

Here we rewrite the standard perturbation power-expansion series for an amplitude in the form of a 



formal integral representation. This representation, in contrast to similar ones in [|9|, ||] and in [14|, does 
not involve integration over the intrinsic momentum k. Even more, it is not related to a Feynman integral 
over the momenta at all. The properties of this representation will be discussed in this section and we 
shall use it in the next section as a convenient "perturbative tool" to interpret the BLM procedure by 
means of the average virtuality flow. 

Let us consider the formal perturbation series s(a s ) for the two-point amplitudes as functions of the 
external momentum Q 2 . The coupling a s is normalized for the default choice of the scale [i 2 = Q 2 . In this 
case, the coefficients of the expansion, d n , are numbers in MS-like schemes due to the cancellation of the 
logarithms ln fc (Q 2 //i 2 ), but they leave behind remnants: the constant parts (ln(C)) of these logarithms 
(accompanied by the /3-function coefficients) are left over in the d n coefficients and have to be taken care 
of. For further convenience, we introduce a new expansion parameter A = boa s "normalized" by the 
factor bo so that 

s(a s ) =d + Y a^d n = S(A) = d + ^ • V A n D n with D x = 1 (2.1) 

and define new coefficients Di = l — - that will simplify the intermediate calculations, and will help 

d\b 

us maintain the contact with the "large &o" limit, bo ^> 1 at A < 1. Note that in the real QCD theory, 
and below the c-quark threshold (for Nf = 3), we have at /i 2 = I GeV 2 the values bo = 9 S> 1 and 

A{n 2 ) = a s ( / u 2 )— — ~ 0.32 < 1 at the NLO level. The running of the coupling A — > A(t) (or a s — > a s (t)) 
follows the renormalization group (RG) equation 

±A = B(A) = -(A 2 + cii 3 + c 2 l 4 + . . .) (2.2) 
th a = 

b, 



with a = tj^j, (2.3) 



o 



where B(A) is the modified /3-function and t = ln(Q 2 /A 2 ) is a natural variable for MS-like schemes. 

At the one loop order we have a well-known solution to ( |2.2| ): A^(t) = — ; at the two- loop order the 

exact solution A^) can be realized in terms of the appropriate branch of the Lambert function W-\, see 
the discussion in Appendix |A], and [17, |l^] for details. Beyond the two-loop level, explicit solutions are 
difficult to obtain. The loop order / is indicated by the subscript I in parenthesis, e.g., A^y 

To evaluate the sum in ( [2.1| ), we need an appropriate representation for the coefficients D n . To this 
end, let us decompose these coefficients into two parts 

D n = D + (n) + (-l)^D-(n) (2.4) 
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to distinguish the values of the coefficient for even and odd values of the index n. Then, we construct 
the generating functions V , V~ , respectively, for these parts D + (n), D~(n) in the following way: 



D ± {n)= / V ± {a)a n - 1 
Jo 



da 



with the normalization condition 

L>i = D + 

At large a the behavior V ± (a) 



poo 

l)+£r(l)= / (V + (a) +V~{a)) da = 1. 
Jo 



(2.5) 



(2.6) 



~ a 



7+1 p-a/c 



corresponds to the expected asymptotic behavior of the 



expansion coefficients D n |L{| at large n 



V ± {a) ~ a 7+1 e 



7+1 -a/c 



D ± (n) ~ D n ~ r(n + 1) n 7 c n , 



(2.7) 



which is a familiar behavior exhibiting a purely renormalon divergence proportional to nl. Below, we 
shall construct a representation for the sum (2T), based on Eq. (^) and on the representation (2.5). 

1-loop integral representation. At 1-loop level, evolution ( |2.2[ ) leads to a useful representation 
for the powers of ^4(i), viz., 



nl 



dt 



A 



(i) 



(2.1 



Substituting ( |2.5| ) and ( |2.8| ) into definition ( |2 . 1[ ) and changing the order of the sum and the integration, 
we obtain the formal integral representation 



S(A (1) ) = do + 



di 
b 



V + (a) 



cxp 



-a- 



d 
"di 



+ 7> _ (a) 



exp a 



d_ 

"~dt 



da 



which contains the shift-argument operators. Finally, we find the representation 

d f°° 

S(A (1) ) = d + j- ■ J {V + {a)A (l) {t-a)+V-{a)A {l) {t + a)}da 
that is linearized in the coupling A\ in terms of the convolution 

(V + (a)A (1) (t -a)+ V~(a)A (1) (t + a)) da = (A {1) (t - a)) + {A (1) (t + a)). 



(2.9) 



(2.10) 



Following Neubert's proposal H, it is convenient to consider the integration in Eqs. (2.9) and ( |2,10D as an 
average of the corresponding coupling A over V . This is denoted as an average of the coupling A(t =p en) 
on the right-hand-side (RHS) of (fTIcD - 

Representation ( p.9j ) seems to be close to that invented by Neubert ||, and also by Beneke, Braun, 
and Ball |], |l0[| , though it turns out to be different. The main difference between the two representations 
will be demonstrated in the next subsection. Note that representation fl2.9j ) is a "formal" one because the 
change of the order of summation and integration in its derivation has not been proved. Nevertheless, 
the integration over the Taylor expansion of the factors At\\ (t — a) and A(\\ (t + a) in the integrand of Eq. 
(p.9j) leads again to the initial series, given by Eq. (2.1), up to any finite order of the expansion. Therefore, 



we shall call this Eq. ( |2.9| ) the "sum" of the perturbative series (2.1) in the one-loop running-coupling 
approximation. The integrand in Eq. (2.9) has a pole singularity at t = a, which reflects the fact that this 
representation is ill-defined. One can now discuss how the contour should be deformed to give a rigorous 
sense to this integral [20, ^]. The residue of the pole can be taken as a measure for the uncertainty of 
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the asymptotic series and it reflects the factorial growth of the perturbation-series coefficients D n (see 

i EL 0)- 

We are not going to discuss here the non-perturbative interpretation of Eq. ([2,9|). Let us only mention 
that quantum field models with an integrable running coupling Ary (t) lead to a finite expression for the 
convolution given by Eq. fl2.10| ). Therefore, the finiteness of the coupling of the analytic perturbation 
theory (APT) |24|], A^ T (t) < 1 provides an example for a convergent perturbation series in the 
sense of the convolution representation (i.e., in a weak sense). For a brief discussion of this feature see 
Appendix |A[ 

i-loop generalization of the integral representation. Let Am be the solution of the RG 
equation in the /-loop approximation on the RHS of Eq. (j2.2j). Then, the initial series, Eq. ( |2.lD , with 
the coupling Am can be represented as 



that should be compared with Eq. fl2,9|). For example, for the two-loop running of A, the representation 



S(A il) ) = d + ^-IJ o V + (a)A {l) U + a-^B-\da+(±^T)\ (2.11) 



given by Eq. (2.10) for the average coupling, should be replaced by 



Representation ( 2.11 ) can be proved in the same way as Eq. ( [2,9| ) in the one- loop case (see Appendix A). 
Its Taylor expansion in the second term of the argument of Am (t =p . . .) generates the PE in Eq. ( |2.1| ) 
for A = Am. In what follows we shall call S(Am) in Eq. ( 2.11[ ) the "sum" of the perturbative series in 



the Z-loop running-coupling approximation. 

3. An Illustration: 1-loop BLM procedure 

Here we consider the standard BLM procedure from the "P-distribution point of view. First, we rewrite 
the series, Eq. fl2.1|) , in the form of the average 

S(A (1) )=do + y-{A (1) (t-a)). (3.1) 

Note here that in the one-loop approximation one can take only one of its parts for simplicity and then 
use the linear property of this average. Expanding the coupling in the average (A^(t — a)) around to 
and taking into account A(to) = Aq, one finds 

(A ((t -t -a)+ t )) =A - A 2 (t -t - (a)) + A 3 ((t - t ) 2 - 2{t - t )(a) + (a 2 )) + .... (3.2) 

In order to fix an appropriate scale to, it is necessary to demand that the NLO coefficient (at the order 
Aq) in the expansion ( |3. 2| ) or some particular part of it (to be defined later) should be nullified. To this 
end, let us consider first the structure of (a) = D 2 in detail. 

D2 is proportional to the coefficient (I2 of the initial series ( pT|) which contains a term proportional 



to bo appearing due to the coupling renormalization. In the discussion to follow, it turns out to be 
convenient to use the following shorthand notation: 

d 2 = d x ■ (b e + f) = d x - ((bo) 1 d 2 [l] + (&o)° d 2 [0}) . (3.3) 
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Here and below, the square brackets contain the powers of the /3-function coefficient bo associated with 
the expansion coefficient d 2 . The coefficient tfcfl] originates from the diagrams shown in Fig. 1(a), (b), 
where the ellipses denote other diagrams with the vertex and the gluon or quark field renormalization. 
The coefficient ^[0] stems mainly from the diagrams displayed in Fig. 1(c), (d), this coefficient is not 
responsible for the coupling renormalization. The notation c^f---] in ( |3.3| ) is particularly useful in the 
multiloop case because it allows the inclusion of the powers of the higher /3-function coefficients. Then, 
we get 

d2 , m , d 2 [0] 



(a) = D 2 



d 2 [l] + 



Two possibilities to fix tg are considered in the literature: 

(i) to = t — D 2 =>■ fj 2 = Q 2 exp(— D 2 ) which means that this coefficient is nullified in Eq. 
content of the coefficient D 2 is included as a whole via the new expansion parameter Aq. 
meaning of the Fastest Apparent Convergence (FAC) procedure || denoted by 



(3.4) 



|). This 
This is the 



(A (1) (t-a)) 



A + A 2 -Q + O(Al). 



(3.5) 



Here the structure of the O(Aq) tail looks like Aq (D3 



n 2 



approximation can be reduced to -Pfac> 



FAC(l) 



+ 0{A{ 
9{±D 2 ) 5(a 



The distribution V in this 



|£> 2 |). 



(ii) Alternatively, we may demand that the contribution d 2 [l] to the coefficient D 2 , which is oc bo, 
is nullified so that to = t — d 2 [l] fJ- 2 = Q 2 exp(— ^[l])- I n this case, only the term ^[1] responsible 
for the coupling renormalization is included by the new expansion parameter Aq. This defines the BLM 
procedure |l], |5| 

(t - a)) = A + A 2 ■ M + {Al). (3.6) 

The remaining term of this procedure proportional to Aq is suppressed by the inverse power of the "large 
6 ". The N 2 LO term in Eq. (Q looks like A^ ^D 3 - 2D 2 d 2 [l] + (^[l]) 2 )- The distribution V in the 

BLM case can be reduced to -Pblm, — > -P^ L m(i) = ^(^^[1]) 5 (a — \d 2 [l]\). This determines an 
intrinsic "scale" ^[1] and shifts the normalization scale from t to to- 
Let us consider now the generating functions V from the point of view of their standard PE. In this 
case, V can be represented as a formal series; viz., 



V ± {a) = ^y^5 {n} (a) • D ± (n + 1) 



(3.7) 



n=0 



that generates just the standard PE, ^2 n= iA n D n , in Eq. (|2.1| ). In practice, only a few terms of this 
series are known. Such a broken series in ( p.ip corresponds to the broken series in ( |3,7j ) where every term 
is strongly concentrated around the origin a = 0. However, in the case of a smooth function V, that is 
not strongly concentrated near the origin, (ET7l) provides only an insufficient approximation to V. Indeed, 







(a) 



(b) 



(c) 



(d) 



Figure 1: Diagrams (a,b,. . .) contribute to b tfc[l]; diagrams (c,d) contribute to c?2[ 



one should take into account much more terms of the expansion in Eq. ([T?]) or even its infinite subseries 
to approximate the real behavior anywhere but not close to the origin, say, near the first extremum of 
the V(a) in a. The first two terms involved into the BLM procedure 

5(a)-5^(a)d 2 [l} + ... 

generate an infinite model sum, like (at d 2 [l] > 0) 

00 (—~\\ n 

V+(a) - P+ LM(1) (a) = 53 H^ W(a) WaM)" = S(a - d 2 [l]) (3.8) 

n=0 

in contrast to the standard PE in Eq. ( |3.7j ) with only a few known terms. The model distribution 
^blm(i) = d(ct — ^[l]) looks preferable because this approximation takes into account at once (admittedly, 
in a rather crude manner) the main feature of behavior of V{a) near the extremum at a ~ d 2 [l]. What 
parts of the PE in the next higher orders should one involve in the procedure to improve Pblm? Or, 
stated differently, which are the diagram classes that generate these parts? To clarify this issue (and to 
define the next approximation of V), one should analyze the structure of the d n coefficients at N n LO of 
the PE. This is the task of the next section. 



4. Perturbation expansion — from series to matrix 

In the previous section, we considered the structure of the coefficient d 2 , Eq. (|3.3| ), which is the basis 
of the standard BLM. Now, let us extend this representation to the next higher coefficient appealing to 
the evident diagrammatic origins of the renormalizations there. In N 2 LO, see examples of the diagram 
elements in Fig. 2(a,b), the renormalization of a s coming from one-gluon exchanges and from vertices, 
and quark field renormalization of the same order of expansion generates contributions proportional to 
o 3 6q, a ?J b\. These terms should be compared with those originating from Fig. 1(a) and Fig. 1(b) in NLO. 
A new kind of contribution like a 3 bo is generated by the two-gluon exchange with the renormalization of 
one of these gluon lines/ vertices, see Fig. 2(c). The representation for d% is similar to that in Eq. ( |3,3j ) 
and looks like an expansion in a power series in bo, b\, . . ., i.e., 

d 3 = df (^ 3 [2,0]+6id 3 [0,l]+M3[l,0]+4[0,0]) , (4.1) 

where the first argument uq of the expansion coefficients d n [no,n\, . . .] corresponds to the power of bo, 
whereas the second one n\ corresponds to the power of 61, etc. The coefficient d n [0,0] represents so- 
called [^] "genuine" corrections with n« = for all possible bi powers. One typical diagram generating 
such contributions is depicted in Fig. 2(e) (compare with the diagrams in Fig. 1(c) and Fig. 1(d)). If all 
the arguments of the coefficient d n [. . . , m, 0, . . . , 0] to the right of the index m are equal to zero, then we 
shall omit these arguments for simplicity and write instead the brief notation d n [. .. ,m]. 

In N 3 LO the renormalization of a s generates contributions proportional to af&o, a s^obi, o^b 2 orig- 
inating from one-gluon exchanges and/or vertices. On the other hand, contributions proportional to 
Og&Q, dgbx, a^bo originate from the mixing of the coupling renormalization stemming from different 
gluon lines and/or vertices. Finally, contributions like appear from "genuine" corrections, i.e., those 
corrections that come from diagrams which do not contain gluon or quark or intrinsic vertex renormaliza- 
tions. Such diagrams do not generate contributions containing powers of the beta-function coefficients. 
Following these lines of argument, the d^ coefficient looks in this brief notation like 

d 4 = d x ■ (6gd 4 [3] + hbodill, 1] + b 2 d 4 [0,0, 1} + b 2 d 4 [2] + M 4 [0, 1] + 60^4 [1] + ^[0]) . (4.2) 
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The same ordering of the /3-function elements holds for all higher coefficients d n . It is convenient for our 
purposes to present this /3-function structure in terms of the "normalized" variables A and D n . The D n 
coefficients have an evident form (in the left column the corresponding coupling powers A n are shown) 



A\t) D 1 

A 2 (t) D 2 

A 3 (t) L> 3 

A 4 (t) D 4 



A N (t) D N = d N [N-l] + 



1; 



1 



d 2 [l} + r -d 2 [0}; 
d 3 [2] +cid 3 [0,l] + 



b 



d 3 [l] + T-d 3 [0] 



d 4 [3] + cid 4 [1,1] + c 2 d A [0,0,1] + 



~r~ ■ (d 4 [2\ + Cl d 4 [0, 1] + ^ ■ (dill] + ^d 4 [0) 
bo \ b V b 

A 5 {t) D 5 = d 5 [A} + c 1 d 5 [2,l]+c 2 1 d 5 [0,2}+c 2 d 5 [l,0,l]+c 3 d 5 [0,0,0,l} + 



1 

Vo 



d5[3] + cid 5 [l,l]+C2d 5 [0,0,l] + 



1 

bo 



d 5 [2] + ci4[0,l] + 



d 5 [l] + -d 5 [0] 



(4.3) 



where Cj are defined in Eq. (^ 

Here, we do not discuss how to derive the elements of this representation explicitly for the known 
multi-loop results. It is a separate task that is solved for the partial case of the Adler function in Appendix 
|D| . Note in this context that similar representations for some physically interesting cases were presented 
in ]i~6|, p5|. In the considerations to follow, we assume that the elements determining the coefficients D n 



in Eq. (4.3) are known so that they can be used for further considerations. To clarify the meaning of the 
structure of the above representation, let us comment on the nature of its elements. The first column 
of the coefficients d n [n — 1] in Eq. (4.3) corresponds to the "bubble approximation" that includes the 



contributions from the diagrams with the maximum number of bubbles in this order, see Fig. 1(a). These 
"bubble" contributions involved in the extended BLM procedure, were considered before in || H [l(j, ||] . 
However, there are other unsuppressed contributions in Eq. ((D|), which have been underlined for a faster 
identification. In fact, for the MS- scheme the values of the coefficients q, namely, 



ci « 0.79; c 2 ~ 0.88; c 3 w 1.9; c 4 = unknown 



(4.4) 





+... 





(a) (b) (c) (d) (e) 

Figure 2: Diagram elements for a s -radiative corrections contributing, e.g., to the photon polarization operator II: (a) 
a chain of ghion bubbles contributing to the &o-term in the coefficient ds. The ellipses here denote other diagrams with 
the gluon/vertex renormalization; (b) diagrams contributing to the &i-term in ds to illustrate the renormalization of gluon 
fields and vertices; (c) diagrams that contribute to the fco - term in d$\ (d) diagrams generating 62- and 6o&i- terms in di\ (e) 
diagrams without any intrinsic renormalization contributions to ds[0]. 
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are of order 1 where, for definiteness, their estimates have been obtained for Nf = 3 so that one has no 
reason to neglect them in Eq. (|4.3|) . These terms (cf. the ci-term in D3 in Eq. ( f4.3|) ) originate in part 
from the diagrams in Fig. 2(b), while the 02- and c\- terms in D4 in Eq. fl4.3| ) stem from Fig. 2(d), and 
so on. To put our final results on a broader basis, one may imagine that the uncalculated coefficients Cj 
are also of the order of unity, q = 0(1), as we discussed in the Introduction. 

We see that the PE series ( |2.1| ) subject to the representation ( |4.3| ) has a complicated structure. We 
face two different expansion parameters in Eq. ([4.3j): first, the coupling A for the lines and, second, b$ 1 
for the columns. To simplify the treatment of these parameters, it is convenient to introduce the notation 
A 1 ■ yij ■ b^ +1 for the contributions in Eq. ( [4.3j ) and Dj = y+j ■ b^ +1 for their coefficients. Y = {yij} 
is a triangular matrix with the elements y nn = d n [0] on the diagonal. These diagonal elements do not 
contribute to the coupling renormalization and are maximally suppressed by the 6q 1 powers in D n , while 
the unsuppressed terms (underlined in Eq. (f4.3|)) are contained in the first column of Y 



yu = I; (4.5) 

2/21 = d2[i]; (4.6) 

2/31 = rf 3 [2] + cid 3 [0,l] ; (4.7) 

2/41 = d 4 [3] + cid4[l,l]+C2d4[0,0,l] ; (4.8) 

All these terms originate from the renormalization of a single coupling/gluon line. As an illustration, see 
the skeleton diagram in Fig. 3(a) which characterizes all diagrams like those in Fig. 2(a),(b),(d). Note 
that only the terms yu survive in the formal "large-60" limit, 60 ^> 1 at c, = O(l). 

5. The BLM task, first stage of the generalization 

Here we construct the first stage of the generalization of the BLM procedure, basing it on the hierarchy 
of the contributions introduced in the previous section. We start with the most important unsuppressed 
contributions yu. The goal is to find a new pair (tx,A(tx)) — following the RG law — which is able to 
nullify all the yu terms and then absorb them into the new expansion parameter A(ti): 

A(t) ^ A(tx) = Ax; ii=t-Ai 

Ai = Ai(Ax) = A 1)0 + Ax • Ai,i + A\ • A lj2 + . . . (5.1) 

The shift Ai from the original scale t to the first intrinsic scale tx is going to be presented in the form of a 
perturbation series in Ax (as first suggested in ]l2|). The corresponding procedure consists in re-expanding 
the initial expression A^yijbQ^ 1 in terms of the new coupling Ax and subsequent rearrangement of the 
power series. This includes the re-expansion of the coupling A(t) = A(Ax, Ax) and its powers in terms 
of A l5 

A(t) = A(Ax,Ax) =Ax- B(Ax)^ + B'{Ax)B{Ax)^ + ... 

= exp(-AxB(A)d A )A\ A=Ai , (5.2) 



which can be realized with the operator exp (— AxB(A)8a) [■ ■ ■] \a=Ai, see, e.g., [22]. Applying this 
expansion together with the expansion for Ai in Eq. ( |5.1| ) to the initial expression A % yijb^ +1 , one 
arrives at the rearranged series in A\, notably, 



e W {-Ax(Ax)B(A)d A ) A*y l3 b^ +1 \ A=A = (Axf $V +1 - ( 5 - 3 ) 
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The elements of the new matrix y^ are illustrated below in terms of the initial coefficients D{ and read 

(5.4) 

; (5.5) 

• D 2 - A ll0 ci + A? i0 - A lii; (5.6) 

• D 3 + (3Af >0 - 2ciA 1>0 ) D 2 - c 2 Ai i0 + ~ciA? )0 - A? + 

(2Ai, - 2D 2 - a) Aa,i - A li2 ; (5.7) 
A N -D N ^ (Axf ■ D N -(N- l)Ai, • Av-i + . . . 

The sequential BLM procedure requires that the y$ contributions should cancel at each order A\ and 



A 


• Di - 


{ A 

* [Ai 


,1 
) 


i; 




A 2 


■D 2 - 


♦ (Ai 


) 2 




- 1Ai, 


A 3 


■D 3 - 


* (Ai 


) 3 


£3 


- 2Ai, 


A 4 


■D 4 - 




) 4 


D 4 


- 3Ai, 



i > 1 in the set of Eqs. (5.5), ( |5.6| ), (5/7),... . This requirement completely determines the elements A^j of 
the expansion of Ai in the original set of algebraic equations. Explicit expressions for the first few Ai k, 
and the omitted term in the equations below, are presented in Appendix |B[ At this point, we compile a 
few coefficients that are important for our discussion: 

(5.8) 
(5.9) 
(5.10) 

Note that the NLO BLM correction in Eq. (|5.9| ) cancels in the particular case d 3 [2] = (d 2 [l]) 2 , d 3 [0, 1] = 
^2(1]. The first of these terms corresponds to the geometric progression of the leading logarithms d^ LM [n — 
1] = (c?2[l]) n_1 , whereas the second one corresponds to the cancellation of sub-leading logarithms. If 
one applies these conditions to Ai 2 in Eq. ( |5,10 ) (see also Eq. ( |B.4| )), one obtains again the evident 



Ai, 


= 2/21 








= 2/31 


- (2/21) 2 - 


C12/21 = dg [2] - d 2 2 [l] + ci (d 3 [0, 1] - d 2 [l]) ; 


Al, 2 


= 2/41 


- 3y 3 i2/2i 


- 2(2/2i) 3 - ci • . . . = d 4 [3] - 3d 2 [l]d 3 [2] + 2(d 2 [l]) 3 + c x 



"geometric" condition ^[3] = (d 2 [l]) 3 for the cancellation of the leading logarithmic part at the next 
step and so on. This test demonstrates the self-consistency of the calculations. If the coefficients d n+ i[n] 
in y n i follow the RG law for the leading logarithms, then they are taken into account already by the first 
element A^o, while in order to take into account the sub-leading contributions, one needs more subtle 
conditions. The absorbtion of such contributions into the proper scales, see, e.g., the c\ -terms in Eq. 
(|5.9|), marks the differences between our procedure and other BLM extensions, for example, that of Ref. 

In this way, one can rearrange step by step for any fixed order N of the PE the first column yn in 
terms of Ai. As a result of this procedure, the initial series can be reduced to a new one that contains 
only one unsuppressed term A\ ■ 1, while all the others are suppressed by powers of 6q , i.e., 

£ (Ay yij ■ b ^ ^ £ (Aif y W • V +1 = A, + Al ■ £ {Aif~ l yff ■ (5.11) 

i>j>l «>j>l «>j>2 

In other words, the matrix Y transforms into a new matrix Y^ 1 ', the first column of which is now 



y^l = Si 1 and the other few elements are presented in Appendix B, Eq. (B.l). These changes in the 
matrix Y (or in the coefficients Di) are absorbed into the new normalization "scale" t\ = t — Ai of the 
coupling A\ = A(t±). Let us stress that no approximation has been used in all steps of the procedure 
up to the order. The above mentioned detailed hierarchy of the contributions is used only in order 
to establish the order of the sequential steps, rather than for performing an approximation. Indeed, 
we initiated this rearrangement with the first unsuppressed yn column, the most important part of the 
hierarchy. Now we will continue with the second column of Y^\ y\ 2 \ that is suppressed by the power 
b^ 1 and proceed in similar manner with this column. 
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6. Sequential BLM procedure, next stages 



Let us continue with the diagonalization of the matrix by rearranging its second column. First, we 
single out the 6q 1 - suppressed terms on the RHS of Eq. ( 5.11[) which are accumulated by its {Y^')r n _i)- 
minor part in 



i>j>2 



(6.1) 



The elements of this minor are represented in the right part of Table |l| (put on the right of the double 
vertical line). The first power 6q — suppressed terms are the first entries in each row of this column. 
Repeating the same procedure as at the first BLM stage (see the previous section) with the column y^ , 



Table 1: The structure of A x ■ y^b^* 1 , here y$ = y^/V: 



22 



Ay 


1 




minor elements 


Ay 


0+ 


(^l) 1 • 


d 2 [0] 
bo 


• 1 






A v 


0+ 




d 2 [0] 
bo 


-(1) ^ <h[0] 

%2 + h2 
°0 






Ay 


0+ 




d 2 [0] 
bo 


' $42 + 12^43'' + 
°0 


<k[0] 
h 3 




Ay 


0+ 
0+ 




d 2 [0] 
bo 


■s£' + i U ) + 

°0 


-v {1) + 

1,3^54 + 
°0 


ds[0] 
bi 



we rearrange again these terms with the help of the new expansion parameter A(t 2 ) at the new scale t 2 
to obtain 



A{h) ^ A(t 2 ) = A 2 ; 



t x - 1 2 = A 2 = A 2 , + A 2 ■ A 2 ,i + A\ • A 2i2 



(6.2) 



like in Eqs. fl5.1|), (^). As a result, all the terms from y^\ except the first one, y^ = £^[0], are 
transferred into the scale t 2 of the new coupling A 2 — therefore, y^ —>■ y^ = y 22 &i 2 - Picking out this 
diagonal term on the RHS of ( |6.1[ ), one arrives at the result of the second stage of the procedure 



(6.3) 



*>i>2 



i>j>3 



To fix the new scale t 2 from A 2im we use the same procedure as in the previous section. The first equalities 
in Eqs. (|5. §D~ fl5.1C| ) remain valid also for the elements A 2>m originating from with an evident shifting 
of all its indices by 1. Therefore, using Eqs. ( [B.lD to determine and choosing the common factor 
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-|— - to normalize the elements (y^ = Vni/V 22 )' one arr i ves a t 



a ~(l) V32 d 3 [l] 

A 2,o = V 32 = 2y 2i = — — - - 2d 2 [lJ; (6.4) 



^2,1 



V22 d 2 [0] 

y$-(y$) 2 -ay$; (6-5) 



Then substituting Eq. (|6. 3|) into the RHS of Eq. (5.11) one obtains the final result of the second stage 



E ( A Y Vii ■ b o J+1 — • • • — M 1 + A % ^ + £ {A 2 r l v$tf +1 ] ] . (6.6) 

i>j>l \ V ° «>i>3 



In this way, one can apply a chain of transformations Y — > — > Y^ 2 ) Y^ -1 ) from column to 

column to complete the diagonalization of the matrix Y, as it is displayed for the first stage in the left 
column of Table |]. At each stage of this diagonalization procedure, one will obtain a new coupling A(ti). 
The final result of this successive BLM procedure reduces the initial perturbation series, Eq. (|2.1|), to a 
series of the special form 1 

S(A) - S'seBLM = do + ^ • d x (l + (l + ^cfc[0] (1 + ...)) j (6-7) 

that contains only the "genuine" elements di[0] = ya, accompanied by the corresponding coupling at its 
proper scale tj, ti = t — Ai — ... — Aj. In other words, the seBLM procedure transforms the standard 
power series a™(t) to the series of the products YYl = ia s (ti) keeping invariant the diagonal elements ya. 
These elements turn out to be the coefficients of the new series expansion. 



Note that Eq. (6.7) can be easily presented in the form 



-S'seBLM = do H 1 1 | ) L-| , (6.8) 

1 a 2 d 2 [OJ 

"l + a 2 d 2 [0] ^ 3[ ° ] 



where aj = a s (tj), a form, motivated by the continued-fraction representation. The benefit of this repre- 
sentation is that it effects at every step the new coupling at its new appropriate scale. This equation itself 
may become the source of a new approximation for the initial series S(A). An instructive interpretation 
of the procedure from the "distribution" point of view is presented in Appendix [C]. The diagrammatic 
sources for both the new expansion coefficients di[0] and the new proper scales t, are discussed there. 
In particular, the skeleton diagram with a single dressed-gluon connection (or vertices) in Fig. 3(a) 
generates the term A(ti)/bo ■ d\ in Eq. ( |6.7D , the shift of the scale Ai is determined by the coupling 
renormalizations appearing from all corresponding partial diagrams, see Figs. 1, 2(a), and 2(b), and so 
on. Representation (|6.7D corresponds to an expansion over the classes of skeleton diagrams and provides 
one with a particular model for V + (a). 



The final result of the seBLM procedure is represented by Eq. (6/7) or Eq. (|6.8|). Let us reiterate 



that the underlying mechanism for obtaining this result is due to the rearrangement of the PE via the 



Let us stress here that the form of this series differs from those suggested in M, 
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y-matrix, which also fixes the hierarchy of these coupling renormalization contributions. The constructed 
procedure takes into account all sources of the coupling renormalization appearing in the PE and absorbs 
them completely into the coupling scales, just in the spirit of the original BLM method Q. At first glance, 
this looks as an improvement but it is rather formal because the seBLM procedure may result in a shifting 
of the scales into a region where the applicability of pQCD may not be justified for both (A(ti), tj) and 
the new expansion coefficients dj[0]. The question is under which conditions this procedure will provide 
us with a better convergent series. The standard BLM procedure may improve the NLO approximation, 
see Eq. fl3.6|), if both coefficients d 2 [l] an d d 2 [0] have the same sign. Similarly, to the above the seBLM 
procedure may improve the convergence of the series when the elements of every line of the matrix Y 
have the same sign. In this case, the absolute magnitude of the final expansion coefficients di[0] (see 
(|6.7l)) is less than that of the initial ones D{. One, however, should not expect an improvement for the 
case when the elements in the lines of Y are alternating in sign and a strong cancellation appears in the 
coefficients Di = y^j ■ b^ +l between different terms. We shall apply the seBLM procedure to the known 
N 3 LO calculation of the D-function in Section ^ and argue that this is an example just for the latter 
case. 



7. seBLM procedure for the D function 

Here we consider step by step the results of the seBLM procedure to highlight its advantages and disad- 
vantages from the point of view of the convergence of the PE of the Adler D-function. 

The initial well-known series for D |2(| can be rewritten by means of the /^-function coefficients to 

read 

D = 3 Yl Q}i d o + di [a + a 2 d 2 + a 3 d 3 + a 4 (f 4 + •••]}, do = 1; di = 3C F ; 
f 

d 2 = b ■ d 2 [l] + d 2 [0]; 

d 3 = b 2 - dg [2] + h • d 3 [0, 1] + b ■ d 3 [l] + ds[0]; (7.1) 
d A = bl ■ d 4 [3] + bob! ■ cU[l, l] + b 2 - d 4 [0, 0, 1] + bl ■ d A [2] + . . . 

However, to recast D into this form is a separate problem which is solved in Appendix ^ on the basis 
of the results obtained in [^| for the MSSM with light gluinos. This new degree of freedom allows us 
to distinguish the contributions with different coefficients 6, in d 3 . Note that the expressions for the 
expansion elements in (7.1) remain valid even including light gluinos that contribute to the /^-function 
(Appendix [A]). The explicit expressions for the elements d 3 [m^n\ are presented in Appendix [D|, while 
below their numerical values are displayed and used 

d 2 = b ■ 0.69 + (7.2) 



d 3 = &§• 3.104 - h-1.2 +b - 55.70 - [573.96 + 19.83 ^ ^ ) . (7.3) 

For the sake of illustration, we substitute the value bo(Nj = 3) = 9, b±(Nf = 3) = 64 into ( |7.3| ), 

d 3 = 251.1 - 76.8 + 501.3 - (573.96 + 0) « 10L9, (7.4) 

in order to compare the contributions stemming from different sources. Furthermore, we shall apply the 
seBLM procedure to D step by step to remove, respectively, the 6o-contribution to N 2 LO, the 6q an d 
^-contributions to N 3 LO, etc., and the associated results will be analyzed. 
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Now we start with our procedure with the original BLM scale setting by transforming the coefficients 
d 2 , d 3 (cf. expressions (|7.2[)-([7^)) and the coupling as follows 



d 2 -h. d 2 = b -0 + ^; (7.5) 
d 3 -> d 3 = bl [d 3 [2] + d 3 [0, l]ci - d|[l] - d 2 [l]c 1 ) + b (d 3 [l] - 2d 2 [0]d 2 [l}) + d 3 [0] (7.6) 

= &§( 2.1555 - 1.0251 ) + &o( 55.70 - 0.46 ) + ■■■« 14.7 (7.7) 

A(t) -» i - ti = Ai, = ~ 0.69. (7.8) 

We see that at the first step (employing the same condition as when we did with Eq. (|^J|)) the value of 
2/31 i s reduced by approximately a factor of 2, whereas the value of &o 2/32 practically does not change. 
On the other hand, the value of the total coefficient is reduced to 14.7 as compared with the initial value 
d 3 ~ 101.9 in Eq. ( |7.4j ). This strong cancellation appears due to the large and negative value of the 
"genuine" term ^[0]. The contents of d^ in Eq. (pM|) also transform following Eq. (5.10). Appealing to 



the results provided in [28], which lead to d^[3] ~ 2.18, one can predict the modification of the "bubble 
part" di[3] of the d&, 

d 4 [3] « 2.18 -> d 4 [3] - d 3 [2]d 2 [l] - 2d 2 [l]{d 3 [2] - d 2 [l} 2 ) » -3.3, 

which looks rather moderate. 

As the next step of seBLM, the modified y 3 ± term in Eq. ( |7.6| ) is transferred into Ai, following Eqs. 



da -> & = &o • + h ■ + b ■ {d 3 [l} - 2d 2 [0]d 2 [l]) + d 3 [0] « -77; (7.9) 
A(t) A(ti); t-t 1 = A 1 = d 2 [l] + A(h) ■ [2] + d 3 [0, l]ci - d 2 2 [l] - d 2 [l]c x ) ; (7.10) 

Ai « 0.69 + A(ti) • 1.13. (7.11) 



Notice that one can put t\ ~ t — d 2 [l] for the A argument in Eq. (7.11) rather than solve Eq. (7.10) with 



respect to t\. The new absolute value of d 3 , d 3 — > ^3 Rj —77 is significantly larger than the value of this 
coefficient at the first step d 3 ~ 14/7. At the same time, the first perturbation correction to Ai in Eq. 
( [7.11 ) looks admissible due to the moderate size of this term. Nevertheless, one can conclude that the 



next step of seBLM does not really improve the convergence of the perturbation series due to the large 
absolute value of d 3 . 

Finally, let us consider the results of the second stage of the seBLM procedure. Following Eq. (|6.4| ). 
one finds that 

1. t\—t 2 = A2,o = d3[l]/d2[0] — 2d 2 [l] « 166 (!); therefore, t 2 t\ — 166 is outside the pQCD domain; 

2. this scale-fixing procedure does not lead to a decrease of the d 3 term d 3 —* d 3 = d 3 [0] m —574 due 
to a large value of the "genuine" term, compared to the contributions of the other terms in Eq. 



Here we encounter the particularly bad case when the second stage of the seBLM procedure cannot be 
applied even to the standard pQCD domain (see the first item above). But, if one, nevertheless, tries 
to apply it, this will not improve the convergence of the perturbation series (the second item) due to 
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I c?3 [0] | S> d 3 . This happens because the procedure ignores the specific properties of the considered series, 
like alternate signs, etc. Let us recall that a similar effect at the LO was mentioned even in the pioneer 
work ||] in considering the branching of the T-decay, T(T — > hadrons)/T(Y — > n + ^). The next section 
is devoted to the solution of this problem. 



8. Generalized BLM procedure to improve the series convergence 



How to improve the seBLM procedure. With hindsight, we now would be tempted to think that it 
would be better to have not performed the second stage at all and try instead another way to optimize 
the value of d 3 after the first one. Indeed, it is not mandatory to absorb the elements yij as a whole 
into the new scale following BLM, but take instead only that part which is appropriate for the coupling 
renormalization. 2 In our case it is tempting not to remove the contribution y 3 i = d 3 [2] + cic^O, 1] 
completely, as we did in Eq. ([Fl]) at the second step above, but rearrange a part of it and absorb one 
part into the coupling renormalization while keeping the other positive part in the remaining expression 
in order to compensate the large and negative d 3 [0] contribution. It turns out that it is convenient to 
introduce the x-portion of 2/31, x 1/31, in order to absorb it into the scale, see Eq. fl3.2| ), while its (1 — x)- 
portion, (1 — x) 2/31, is kept explicit in order to be able to cancel the negative "genuine" d 3 [ti\ in (|8.l|) . 
This trick leads to a modified seBLM procedure termed x-dependent BLM, (xBLM), 



d~3 -> da = &§■(!- x)(d 3 [2] + Cl d 3 [0, 1]) + b ■ {d 3 [l\ - 2d 2 [0]d 2 [l}) + d 3 [0], 
- tl = A a = d 2 [l] + A(h) ■ (x(d 3 [2] + d 3 [0, l]ci) - d 2 2 [l] - <fc [!]<*) . 



(8.1) 
(8.2) 



Let us now establish an "optimization" condition, similar to Eq. ([3.5|), i.e., put d 3 = to determine the 
value of x. One has to make sure that in the result the perturbative corrections are improved for both 
d 3 and Ai (see the fifth and sixth columns in Table ^) in comparison with those ones in Eqs. (^|) and 
( [7.11 ). The final result for D is then reduced to 



D = 3j2Q}i 1 + 3C F 

/ 



a(ii) + --a / (i 1 ) + 0-a d (ti) 



i.3) 



As an illustration, the estimates of A\(Qf) at Q\ = 3 GeV 2 ; Q\ = 26 GeV 2 are shown in the sixth 
column of Table |2| 



Table 2: Results of xBLM for d 3 = ( f 3 = ) and Q\ (sj) = 3 GeV 2 ; Q\ (s 2 ) = 26 GeV 



N f 


60 




X 


t - 


h = Ai 


Ai(Q 2 ) 




X 


t(s) 


-t(h) = 


Ai 


Ai(s) 


3 


9 




0.56 


d 2 [l} + 


a(h)b ■ 0.18 






1.84 


d 2 [l] 


- a(h)bo 


•3.1 




4 


25 

y 




0.24 


d 2 [l]- 


o(ti)6o -0.45 


Ai(Q?) 
ps 0.58, 


f 3 = 


2.56 


<k[l] 


- a(si)bo 


•3.7 


Ai(ai) 
ps -0.27 


5 


23 

y 




-0.11 


d 2 [l] - 


o(ti)6o • 1-19 


Ai(Qi) 
ps 0.52 




3.63 


d 2 [l] 


- a(si)b 


4.48 


Ai(a 2 ) 
ps -0.05 



This possibility has already been used for the LO BLM procedure, applied to the pion form factor in J29|] and termed 
there BLM, see p. 13. 
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It is instructive to apply a similar procedure also to the observable quantity R(s) — 



associated with D: 



cr(e+e~ — * 



_2 



R(s) = D{s) - d x — ■ b 2 a 3 = 3 V { 1 + 3C F [a + r 2 a 2 + r 3 o 3 ] } , (8.4) 
6 f 

IT 2 

where n = d\, r 2 = d 2 , ^3 = ^3 • fe 2 , see > e -S-> P^C ? ^( s ) = hi(s/A 2 ) and a = a(t(s)). The large and 

3 

negative 7r 2 -term, arising due to the analytic continuation, makes r% also negative (cf. with Eq. (7J3)). As 
a result of the xBLM procedure, one should replace r 2 by <i 2 and by 7*3, whereas the x-dependent term 
in r% transforms to 6g • (1 — x)(d^[2] + cids[0, 1] — vr 2 /3). Then, to obtain a positive compensating term to 
cancel d^ [0] at the next step of the transformation (that leads to f 3 = 0) , one should take x > 1 , consult 
the eighth and ninth columns in Table |2[ The result of this procedure is exemplified by showing partial 
values of shifts Ai(si j2 ) that are placed in the last column of this table. For these examples, Ai < and, 
therefore the optimal scale s\ becomes larger than the initial one s. 

Having developed the xBLM procedure, we are in the position to gain control over the size of the 
coefficients dj +2 and the magnitude of the associated scales tj to balance their contributions. This possi- 
bility is of paramount importance because we need to control several different perturbation expansions in 
the generalized BLM procedure: one is the size of the expansion coefficients dj +2 like ^3 in Eq. flS.l|) , the 



others are for the scale shifts Aj like Ai in Eq. (|8.2| ). Finally, in the result of this procedure, Eqs. 
fljQl) , the coefficients of a 2 reduce to d 2 [0] = - (as for the standard BLM procedure), the coefficients 

(^3) of a 3 , reduce to 0, while the values of the new normalization scale t\ (t(§i)) are moderate and 
under control. This effect - to hold the control over both the types of expansions - is the main advantage 
of the rrBLM procedure, while the condition d% = used here is not a mandatory requirement for it. 
Instead of that one can require to significantly reduce the size of 0(3 for better convergence. This trick to 
improve the convergence of the truncated series can be further generalized. 

Matrix-form generalized BLM procedure. To generalize the seBLM procedure in the way 
mentioned above, let us introduce a lower triangular matrix X = {xij}, Xa = associated with the 
matrix Y instead of employing a single scalar parameter x. The element x^ determines those portions 
of the contribution yij which should be absorbed into the coupling renormalization, according to the 
operation x%j yij. On the other hand, the remainder of the contribution x%j (where Xij = 1 — x^) is 
kept in the PE-coefficient D{. The choice X = brings us back to the initial series before performing 
any transformations, while X = {x 2 i = 1, Xj>j = 0} corresponds to the standard BLM procedure. 
The matrix X = {xj>j = 1} corresponds to the seBLM procedure, whereas all off-diagonal elements 
yij should be absorbed in new couplings. For the latter case, the first column xn = 1 itself leads to 
the first stage of seBLM, as one infers by comparing the second column in Table ^ with Eq. (|B.l[) . For 
the discussed xBLM procedure the first column transforms into {1, Aiy2ix%i, A^y^ix^, . . .} instead of 
the diagonalized structure {1, 0, 0, . . .} in seBLM. The schematic sketch of the first stage of xBLM 
is demonstrated in Table ||| and Eqs. ( |8.5| )-( |8?7|) and one should compare them in conjunction with the 
seBLM procedure, presented in Table [l] (or in Eq. ( |B. 1| ) ) . 

A(t) -» A(h) = Ar, t - t x = Ax = Ai, (X) + A x ■ A 1A (X) + .... (8.5) 

The formulae for the Aj proper scale, Eqs. (|5.8| ), ( |5.9|) ,... remain valid in this case by taking into account 
the obvious changes yn — > yuxn, 



A 1A (X) = y 3 ix 3 i - 2(y 2 ix 2 i)y 2 i + (y 2 ix 2 i) 2 - (y 2 ix 2 i)ci. 



(8.6) 
(8.7) 



Table 3: The structure of A 1 ■A[~ 1 y^\x)b Q j+1 at the first stage of the xBLM procedure 
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1 




Ay 


(^1 


I 1 2/21^21 


+ A l 


2/22 

bo 






Ay 


(^1 


) 2 2/31^31 


+ {A, 


, 2 2/22 


— - 2(y 2 i X21) ) 

V2/22 J 


,(A \2 2/33 


Ay 


+ 

(^1 




+ ... 
+ ... 






+ ... 
+ ... 



Additional free parameters x^ give rise to a total amount of parameters n(n — l)/2 in N n LO of the PE. 
This allows one to perform a "fine tuning" of the coefficients of the series minimizing those expansion 
coefficients pertaining to the considered order of the PE. A more complicated structure of the final PE 
series is the price one has to pay for such an improvement of the convergence of the series. 

At the second stage of the xBLM procedure the contribution from the second column y, 2 transforms 
into the sum ^(^2) 1-1 (2/12^2) /&o that should be compared with the single term (A2) l ~ 1 5i2 2/22/^0 in the 
seBLM one. Here A2 is fixed by the modified condition for A2 

A{h) -> A(t 2 ) = A 2 ; h-t 2 = A 2 = A 2 ,o(X) + A 2 ■ A 2 1 (X) + ..., 



-t 2 = A 2 



2/32 0/ \ 

2^32 - 2(?/21 3?2lJ 

2/22 



with still higher terms being omitted for simplicity and whereas the whole expression is similar to the one 
in Eq. ( |6.4[ ). All other portions yi2 Xi2 of this column should be absorbed into the new coupling A 2 by 
means of the expansion coefficients A 2 j. Further stages of the implementation of the xBLM procedure 
are similar to the first ones and result in 



S(A) - S xBLM = d + ^ • d! 

»0 



1 + (A( tj )) 

i>3>2 



i-1 



x ij Vij ' °0 



J.9) 



In Eq. ( p.9[ ) we lose the clear diagonal form, like Eq. ( |6.7D , but obtain more flexibility to fix the values 
of the A l Di terms in the expansion. The elements of the X matrix for various BLM procedures, are 
collected in Table g. 

Table 4: The stages of the generalization of the BLM procedure 



"fine tuning" matrix 

x^ 


Standard PT 


BLM 


seBLM (first step) 


seBLM 


xBLM 





fe^ij 


6{i > 1) 5u5 X j 


0{i> j) SuSjj 


9(i > j) 


results for PT series 


Eq. Q 


Eq. O 


Eqs. O, (P) 


Eqs. Q, @ 


Table |; 
Eq. Q 



Let us pause for a moment to make some clarifying remarks to a few partial cases: (i) The FAC 



setting in N 2 LO corresponds to the condition (Ai) 2 ( 1/21X21 + 



2/22 

bo 



which has been considered in 



Eq. (3.5). (ii) On the other hand, if we restrict ourselves, say, to the N 3 LO, then we have a matrix X 
with 3 parameters, x 2 i,xsx,X3 2 to optimize the contributions of A 2 D 2 and A 3 D 3 , respectively, 



A 2 D 2 — > A x 



A 3 D 3 ^ A x 



A - A 2/22 

^1^21^21 + A 2 — 

bo 



C 2 ; (8.10) 



{A x ) 2/31S31 + (A2) -r- + {A 2 ) Tr 

bo V V22 J &o . 



C 3 ; (8.11) 



where Ai = A(t - Ai), A 2 = A(t - Ai - A 2 ), (8.12) 

and C2, C3 are appropriate for us numbers, e.g., C2 = C3 = 0. The case discussed in Sec. || in conjunction 
with Table |2| corresponds to the partial solution of the above equations for C3 = 0, A\ = A 2 with x 2 \ = 1, 



X31 = x, £32 = 0. A complete set of solutions to Eqs. ( |8. 1Q|) ( |8.12|) for C 2 = C3 = with respect to 
can be obtained and analyzed numerically. 

9. Conclusions 

In this paper we have considered the generalization of the BLM procedure with the aim to (i) extend 
it sequentially to any fixed order of the perturbative expansion in QCD and investigate the property of 
this extension, (ii) improve this generalization with respect to a better convergence of the perturbative 
expansion. To this end, a new hierarchy of the contributions to the expansion coefficients of a two- 
point function was designed which exploits the following properties of the beta-function coefficients: 
hi ~ 6q +1 . Employing this new hierarchy, we constructed a sequential and unambiguous generalization of 
the BLM procedure [|| to any fixed order of the perturbative expansion beyond the so-called "large-60" 
approximation, a procedure we called seBLM. This seBLM procedure leads to a new expansion series in 



terms of a set of new couplings cij = a(tj), Eq. ( |6.7l ), depending on a set of proper scales tj that can also 
be rewritten in a form close to a continued- fraction representation, given by Eq. (|6.8;). The advantages 
and disadvantages of the seBLM procedure were exemplarily discussed in terms of the four-loop Adler 
D-function. We found that the seBLM procedure may still fail in improving the perturbation-series 
expansion for the D-function at its second stage, just for the same conditions under which also the 
standard BLM procedure fails. To improve the convergence of this series, the seBLM procedure was 
equipped by additional free parameters, encoded in a parameter x (see Eqs. (|8.1|)-( |8~3D in Sec. |8|) in the 
spirit of the Fast Apparent Convergence procedure ||]. These free parameters enable us to reduce the 
expansion coefficients of higher orders in a systematic way and at the same time to control the value of 
the couplings a{tj). This improved procedure has been successfully applied to the series of the D- and 
R(e+ e -^h) _ functions with results presented in Table The final generalization of the BLM procedure 
is provided by the xBLM one, which is now parameterized by a whole matrix X = {x^} of additional 
parameters in order to achieve a better convergence of the perturbative expansion, as it was briefly 
discussed in the second part of Sec. ||. This way, one is enabled to control both the size of the expansion 
coefficients and the magnitude of the coupling aj — and also of the associated scales tj — from the most 
general point of view. 
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A. Proof of the convolution representation 



1, Let A = An\ be a solution of the RG equation with the /-loop B-function. Let us introduce the 

A 2 d - 
operator D t = — - s — . Then, we find D+A = A 2 and 
P B(A) dt 

^(b t jA = (Ar\ (A.i) 

Substituting ( |A.1| ) and Eq. fl2.5| ) into Eq. ( |2.1| ), and collecting the corresponding terms in the exponent 
for, e.g., D + (n) part of the series, 



roo I I ttDt I \ /-co , . 

Y J D + {n)A n = V + {a) ~ r~ A ^ da = / [ exp \oiDt)A(t) 

n=l ^° \m=0 m " I ^° 



da 







V+(a) Ai [ t + a- — )da, (A.2) 



one arrives at expression ( |2.11| ) . In the two- loop case the exact solution can be expressed in terms of the 
Lambert function W(z), |3(], [L7| defined by 

z = W{z)exp(W(z)) . (A.3) 

This solution has the explicit form 



where z t = (1/ci) exp [— t/c\ — 1 + in] and the branches of the multivalued function W are denoted by 

Wk-, k = 0, ±1, A review of the properties of this special function can be found in pQ ]; see also Jl7| . 

2. In the same way, one can obtain this representation for the non-power expansion series {j4 apt } 
that appears instead of the standard powers {A n } in the analytic perturbation theory (APT) |23|, p4| ]. 
The APT expansion of S(A) in ( |2.1| ) can be written as 

5(A) = d + £ • V D n A- -> S APT = do + £ • V D n A APT (A.5) 

The one-loop APT coupling constant is a bounded function of t, A APT (f) = — — < 1. The solutions 

for higher values of the index n of the one-loop constants A APT are not powers of A APT , but they can 
be obtained from A APT by means of the same relation expressed in Eq. (|2.8|), see [24|, namely, 



A^(t) = 1 ("|) n ^ PT (0 + K PT W)" +1 . (A.6) 

This leads to the same convolution structure given by Eq. (2.9) with j4 APT (t) entering into the integrand. 
The finiteness of A APT (i) guarantees the convergence of this integral representation. For this reason the 
corresponding series in APT, represented by Eq. ( |A.5[ ), are really well-convergent [[H]]. 
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3. The required /3-function coefficients with the Minimal Supersymmetric Model (MSSM) light 
gluinos [^] are, 

b (JV>, N g ) = H C A - 1 ^TftNf + ; (A.7) 

61 = y^i - f c A (r R N f + - 4 (W^ F + ; (A.s) 

b 2 (N f , N g ) = 2 ^C\ - NfTji (^C\ + 2 f CA C F - 2Cl) + (N f T R f (^C F + ^C?) - 

9 -§N g C A (Cl) + N g C A N f T R (jC A C F + ^Ci) + (N g C A )^C A . (A.9) 

The 63 coefficient, which includes the MSSM light gluinos, is not yet known, so we present it here in the 
standard |33| simplest form 

h(Ns) . (™5 + s^ 3 ) - (= + - C3 ) Nf + (i~ + - 6 ) „ + 

^ (A.0) 



B. Explicit formulae for and A i fc 

The explicit expressions for the few first elements of the matrix are 



2/22 

bo'' 

1 / n . 2/33 

7- (2/32 - 2y 2 i 2/22) + to-; 

00 »o 



A 1 


- 




i; 


A 2 


£> 2 " 


+ A\. 


+ 


A 3 


■D 3 - 


,A\- 


+ 


A 4 


■D 4 - 


,A\- 


+ 



&0 



1 J/44 
71 (2/43 - 32/21 2/33) + TT' 
% % 



A n -D n ^ Al - + i(y n2 -...).-. • (B.l) 
The explicit expressions for the elements of the proper scales Ai and A2 are given by 



Ai, = 2/2i = d 2 [l]; (B.2) 
Ai,i = 2/31 - (2/21) 2 - c\V2\ = 

= dg [2] - d|[l] + ci (d 3 [0, 1] - d 2 [l]) ; (B.3) 

73 

d 4 [3] -3d 2 [l]d 3 [2] + 2(d 2 [l]) 3 
ci (d 4 [l, 1] - 3d 3 [0, l]d 2 [l] + ^(d 2 [l]) 2 - d 3 [2] j + 

c 2 (d 2 [l] -d 3 [0,1]) + c 2 (d 4 [0,0,1] -d 2 [l]); (B.4) 



^1,2 



3 

2/41 - 3?/3i2/2i - 2(y 2 i) 3 - C12/31 + -ci(y 2 \) 2 + (c 2 - 02)3/21 
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A 2 , = ^-2, 21 = §§-2d 2 [l]; (B.5) 
2/22 d 2 [0\ 

a 2/42 „ 2/32 . c 2 n a 2 a 

A 2 ,i = 3y 2 i h 5y 21 - 2y 3 i - A 2 - ciA 2 , ; 

2/22 2/22 



C. Distribution sense of sBLM. 

It is instructive to illustrate the discussed seBLM procedure from the "distribution" -point of view, which 
is based on the convolution representation provided by Eq. ( |2.11| ). This has been demonstrated in Eq. 
( |3.2| ) for the one-loop running coupling explicitly. One can re-derive the results of the seBLM procedure 
at the first stage using these terms, e.g., for the V + case. This goes as follows. First, one expands 
the factor Am in the integrand in Eq. ( |2.11| ) in the variable a. Second, following the line of argument 
underlying Sec. [5], one can arrive at the same Eqs. ( |5.8| )-( ^T0| ) for the elements. To achieve this goal, 
on should specify a partial generating function Pi describing the distribution of the elements of the first 
column of the matrix Y, J °° Pi(a)a n ~ 1 da = (a n ~ 1 )\ = y n \, and then substitute Pi in the place of V at 
the final stage of the procedure. As a result of the first seBLM stage, i.e., Eq. ( 5.11| ), the convolution 
(A( 2 )(t — ...)) i reduces to 



poo 

A 1 = A(t-A 1 )= / 5 (a - Ai)A(t - a)da 
Jo 



10 

Therefore, the distribution Pi reduces to 5 (a — Ai), an expression resembling the standard BLM result 
6 (a — c?2 [1] ) in Sec. This term Ai incorporates all the "a s -renormalizations" associated with the 
single dressed-gluon propagator with two dressed vertices and the connected quark fields in the skeleton 
diagrams in Fig. 3(a). The common normalizing factor d\ originates from the first tree diagrams of these 
skeleton Feynman graphs. At the second stage of the seBLM procedure we deal with the contribution 
from the second column of Analogously, a normalized distribution P 2 can be introduced for y^ 2 , 

f °° P2(ct)a n ~ 1 da = (a n ~ 1 )2 = y„ 2 . These contributions correspond to "two-gluon" skeleton diagrams 
illustrated in Fig. 3(b) (see also Fig. 2(c)). By now the effective scale of the coupling A(t — Ai) is already 
fixed at the previous stage. As a result of this operation, P 2 — > 5 (a — A 2 ) and the normalization scale 
of the couplings appear to be shifted, A(t — Ai) — > A(t — A± — A 2 ). The normalizing factor d 2 [0]/6o for 
that contribution originates from the undressed diagrams, which correspond to the skeleton graphs in 
Fig. 3(b). 



VA 




(a) ( b ) 

Figure 3: The elements of skeleton diagrams are presented: An oval box denotes a dressed vertex; a thick wavy line 
represents a dressed gluon propagator and a thick fermion line stands for dressed quark propagators; the arrow indicates 
the value of the seBLM "scale" of the corresponding effective charge, (a) Skeleton diagrams for the first stage of seBLM; 
(b) Skeleton diagrams with two dressed gluon lines for the second stage of seBLM. 
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Executing a number of seBLM stages for (Ar 2 \(t — ...)}, one arrives at the representation 
(A (2) (t -...)) da [6(a - Ax) + ^A(t - Ax) • 6(a -A 1 -A 2 ) + ..?j A(t - a). (C.l) 

The kernel of this convolution can be compared with the PE representation in Eq. (|3.7| ), 
(— 1 ) n 

V+(a) = ^-=f-S^{a) ■ D+ ( n ) " A i) + d 2M/b A(t - A x ) ■ <5(a - Ai - A 2 ) + . . . (C.2) 

n 

The final seBLM series, Eq. ( |6.7[ ), corresponds to the expansion encoded in the skeleton diagrams. The 
coefficients of this expansion d n [0] (y nn ) originate from the first undressed skeleton diagrams, whereas 
the scales A n appear due to the renormalizations in these skeleton diagrams. 



D. Representation for the D-function 



The Adler D function is known f27J for the MSSM with N g light gluinos, D(a,N f ,N g ). On the other 
hand, one can obtain explicit expressions for the functions Nf = Nf(bo,b\) and N g = N g (bo,bi) solving 



the set of equations (|A.7| ), (|A.8| ) with respect to Nf, N g . Substituting these solutions into D(a,Nf,N g ), 
one arrives at the expansions Eq. (pi]), (p.2|) and Eq. QD.4[)-(p^), 



D(a) = 3 ^2 Q) { rf o + d\ (a + d 2 a 2 + d 3 a 3 + •••)}, d = 1; d\ = 3C F ; 
f 

d 2 = b d 2 [l)+d 2 [0}; (D.l) 

d 3 = b 2 d 3 [2] + 6i • d 3 [0, 1] + b Q d 3 [l] + da[0]; (D.2) 

d 4 = bld 4 [3) + b 2 ■ d 4 [0, 0, 1] + & M 4 [1, 1] + 6gd 4 [2] + M 4 [l] + M4[l] + <U[0] • (D.3) 



The iV^-terms of cf 4 have recently been calculated in [^8|, but, unfortunately, these results cannot be 
used in our approach. The reason is that it is impossible to separate the terms ^d^O, 0, 1] and bobid^l, 1] 
that are of the order of 0(6q) from the &Q-term, i.e., 6gd 4 [2], that also contributes to the "iV| projection" . 
This is why one needs a new degree of freedom (like the gluino). 

d 2 [l] = y-4C 3 « 0.691772; d 2 [0] =9±-9l. = L (D.4) 

302 76 101 
d 3 [2] = — -— C 3 ~ 3.10345; d 3 [0, 1] = — - 8C 3 « -1.19979; (D.5) 

d 3 [l] = C A Q + yCs - yC 5 ) - C F (18 + 52C 3 - 8OC5) « 55.7005; (D.6) 

1 ^,n2 , ^ , 1 ^2 ) _ 72C 2, , 5 ^176 loo ^ (E/<5/) 2 



d 3 [0] = ^(523Ci + 852C A C F - 4UC F ) - 72(7^3 + — \ — - 128C 3 J 3 g q2) i D.7) 

2 

-573.9607 - 19.8326 ^ . (Di 
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